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ABSTRACT 

Context. GCIRS 3 is the most prominent MIR-source in the central parsec of the Galaxy. NIR spectroscopy has failed to solve the 
enigma of its nature. The properties and peculiarities of extreme individual objects in the central stellar cluster contribute to our 
knowledge of star and dust formation close to a supermassive black hole. 

Aims. We initiated an unprecedented interferometric experiment to understand the nature of GCIRS 3, where we investigate its 
properties as a spectroscopic and interferometric reference star at 10 yum. 

Methods. VLT/VISIR imaging separates a compact source from diffuse, surrounding emission. The VUTI/MIDI instrument was 
used to measure spectroscopically resolved visibility moduli at an angular resolution of ~ 10 mas of that compact 10 ^m source, still 
unresolved by a single VUT. Recent NIR/MIR photometry data were added to enable simple SED- and full radiative transfer-modeling 
of the data. 

Results. The luminosity and size estimates show that IRS 3 is probably a cool carbon star enshrouded by a complex dust distribution. 
Blackbody temperatures were derived. The coinciding interpretation of single telescope and interferometric data confirm dust emission 
from several different spatial scales. The interferometric data resolve the inner rim of dust formation. Despite observed deep silicate 
absorption towards GCIRS 3, we favor a carbon-rich circumstellar dust shell. The silicate absorption most probably takes place in 
the outer diffuse dust, which is mostly ignored by MIDI measurements, but very observable in complementary VUT/VISIR data. This 
indicates physically and chemically distinct conditions of the local dust, changing with the distance to GCIRS 3. 
Conclusions. We have demonstrated that optical long baseline interferometry at infrared wavelengths is an indispensable tool for 
investigating sources at the Galactic center Our findings suggest further studies of the composition of interstellar dust and the shape 
of the 10 fim silicate feature in this extraordinary region. 

Key words. Galaxy: center - AGB: dust shells - Techniques: interferometric 

1. Introduction pergiants ( T < 10^-^ K) prov i de about 50% of the flux ionizing 

, 1 nr:7^7^ . flie reg ion dBlum et a .1119951 iKrabbe et al.lll995t iNajarro eTaLl 

At a distance of -7.6 kpc dEisenhauer et al. | | 2005D the center ggg^ ISchrMer et al.1 (12003 ') have shown that a few tens of 



. of flie Milky Way is by far the closest center of a large spiral carbon-rich supergiants can produce about 50% of the mass loss 
galaxy. Its astrophysical properties can be studied on a unique ^ j^^g^ ^j^jj^ ^^^^pj^ ^jjj^ ^^j^^ neighborhood characteris- 
angular scale of ~ 40 mpc/arcsec, which is two orders of mag- ^-^^ ^j^^^^ ^-^^^^ ^1^^ -^^-^^ ^^^^ ^j^^ ^^^^ j^^^ 
nitude smaller than the angular scale at the nucleus of M31, formation properties of a stellar cluster can be dominated by 
the next comparable galaxy ( iMcConnachie et al. ||2005^. Star for- ^ prominent stars. These facts underline the importance of 
mation and the kinematics of the central stellar cluster can be ^^^^y^^^ individual extreme objects like GCIRS 3 ' to under- 
studied in the region of direct influence of the supermassive ^^^^^ properties of the entire surrounding stellar environment. 
black hole (SMB H) at the dynamic center of the Milky W ay rp, ^ , ^ f j - f A2/\/nn\- , , o 
teckart & Genzellll996tlCxhez et alJl200dlSchodel et alM , ^he recent advent of mid infrared (MR instruments on 8 m 

^ — — : — -r' . : ; . ^ ^ class telescopes enables us to study in detail the thermal dust at 

. I Jf^'^f history of star formation, despite the tidal forces ^j^^ q^j^^^.^ ^^^^^^ ^^^^ unprecedented angular resolution, 

of the SMBH, IS manifested in the existence of numerous mas- investigation of the circum- and interstellar dust distribution 

1 ^^^^' SMS i n the centra l c^ter ( Krabbe et al. 1995; ^^e GC uncovers stellar mass loss, zones of wind interaction, 
iGenzel et alj 120001: lEckart et all 120041: iMoultaka et al.. ,2004.) 



The 7 most luminous (L > 10^'^^ Lq), moderately hot blue su- 



formation history, evolution, and kinematics. 

Photometric and spectral prope rties of dusty stars a t the 

, a- . , TZ Tit n « t\ir\ii ^ r-.u , GC have recentlv been published bv Moultaka et alj (|2004|) and 

Send offprint requests to: J.-U. Pott, now at W.M. Keck Observatory: i , , 'n i^Wr^^ l ArTrTA • ■ — \ — ■ 

jpott@keck.hawaii.edu IViehmann et all ( I2005L 120061) . Despite an average optical extmc- 

* Based on observations obtained at the European Southern 

Observatory, Paranal, Chile (programs 073.B-0249, 075.B-0113, ' in the following IRS 3 

077.B-0028). " 3-20 fim covering the atmospheric L, M, N, and Q-windows 
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Galactic Center at 8.6 /im 




Merging of individual reductions 



Fig. 1. Flux-calibrated Lucy-Richardson deconvolved image at 
8.59 yum after restoration with a 250 mas Gaussian beam. The 
pixel scale is 75 mas per pixel. The logarithmic contours levels 
are 1.6"-7 mJy. The flux scale is given on top. IRS 3 and Sgr 
A* are highlighted. D etails of the data reduction are given in 
ISchodel etal.ll2007bl) . 



tionofAv ~ 25 dScoville et all 120031; IViehmann et al.l l2005h . 
near-infrared spectroscopy and imaging reveal the nature of the 
underlying dust-embedded stars in most cases, since the dust ex- 
tinction is wavelength-dependent and decreases from the optical 
towards longer wavelengths. 

This article focuses on the most prominent of the MIR 
bright dusty sources, IRS 3, the embedde d stellar source, which 
still eludes any spectral classification dPaumard et al.l l2006t 
iTanner et al.l 12006). It is located within the central 20". A re- 
cent state-of-the-art narrow-band MIR image at 8.6 fim with an 
angular resolution of only 250 mas is shown in Fig. [T] The ex- 
tended and diffuse dust emission surrounding IRS 3 is visible at 
unprecedented angular resolution in this image. 

The NIR extinction studies reveal a spatial variation of only 
~10 % (Ak ~ 3) of the interstellar extinction over this re- 
gion (IScoville et al.ll2003t iMoultaka et alj 12004 ISchodel et alj 
l2007ak 

In contrast, narrow-band, A^-band photometry and spec- 
troscopic observations are interpreted to indicate a signifi- 
cant amount of additional 10 fim silicate absorption along 
the line of sight t owards IRS 3 with respect to other parts 
of the central 20" dBecklin et al. 1978; Roche & Aitken 1985; 
IViehmann et al.ll20()6l) . The published intrinsic optical depths of 
about T9. 8(IRS 3) ~ 1, in additi on to the average t9.8(GC) ~ 
3.5 (e.g. iRoche & AitkenI 1 1985b . still underestimate the true 
value due to source confusion problems. While our new high- 
resolution VISIR data clearly indicate that more than 50% 
of the A^-flux is diffuse, extended emission at 0.3" resolu- 
tion, one byproduct of the MIDI observations is the calibrated, 
low-resolution spectra of the compact emission, which show 
T9.8 (IRS 3) to be much larger than the value given above. 

iRieke etaTId 19781) found the spectro-photometric MIR prop- 
erties of IRS 3 to resemble either young stars or OH/IR stars. The 
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Fig. 2. Part of the science data taken on 2005-05-25. We overplot 
the different reduction methods (coherent and incoherent fringe 
averaging and static and dynamic detector apertures) to demon- 
strate the different results. The labelling indicates the artificial 
wavelength shifts successively applied to increase the readabil- 
ity of the plot. The relative errors of each method is given. The 
rightmost error bars present the merged, maximum error bars. 
We based our analysis on these merged data. 



latter interpretation is opposed by the lack of OH-maser emis- 
sion. With an MIR color temperature of ~ 400 K, the central 
emission of IRS 3 was found to be (together with the nearby 
GCIRS 7) the hottest and most compact of th e sources dominat- 
ing the thermal du s t irradiation from the GC dRieke et alJI 19781 ; 
iGezari et al]|l985l; [Smith et al.lll990l) . Extended dust emission 
around IRS 3 interacting with external stellar winds has been 
found in recent L- and M-band observations ( Viehmann et aT| 
200 5J). While hot star hypotheses are given by some authors 
( Krabbe et al. lfT995t iTanner et alJl2003 1 the lack of ionizing gas 
leads Roche & Aifl^cenl dl9'8¥ to the assumption of IRS 3 being 
a cool dust-enshrouded star. 

Within the past two years we collected a unique dataset of 
optical long-baseline interferometric (OLBI) data of IRS 3 at 
10 fim. These constitute the first successful OLBI observation 
of an object within the central parsec of our galaxy, opening 
the window to NIR/MIR GC observations at highest angular 
resolution (Pott et a l. 2005a). In this article we show that the 
OLBI data strongly support the former hypothesis of a cool 
dust-enshrouded star. Furthermore our results shed light on the 
amount and spectral shape of the interstellar 10 yum extinction 
towards IRS 3 at unprecedented angular resolution. Since MIDI 
is a relatively new instrument, the achievable precision of vis- 
ibility measurements of such distant, challenging targets is not 
common knowledge yet. In Sect. |2] we therefore describe the ob- 
servations, the extensive data-reduction and calibration efforts, 
and the evaluation of different data reduction techniques in some 
detail to ease and stimulate similar experiments and the compar- 
ison of their outcomes. Then the immediate observational results 
are given (Sect. O, followed by a detailed discussion of the re- 
sults in the astrophysical context (Sect.|4]i and a summary of our 
conclusions in Sect|5] 



2. Observations and data reduction 

In 2004 we started an observing campaign to study the bright- 
est MIR-excess sources in the central parsec with MIDI at the 
highest angular resolution available today. The MID-Infrared 
interferometric instrument (MIDI) combines the light of two 
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Table 1. Observing log of IRS 3. The wavelength dependent TF 
has been estimated for each night on the basis of regularly con- 
ducted calibrator measurements, typically about once per hour 
(Sect. |2]i. The number of calibrators used per night is given, and 
the applied parameters are given in Table |2] 
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Baseline 


PB 


PA Airmass Accuracy 




[m] [degEofN] [1] 


[%] 




#1 night: 2004-07-07 


5 calibrators 






U2-U3 


45.7 


46 1.03 


12 




U2-U3 


42.3 


J A, i . i U 


12 




U2-U3 


37.0 


55 1.45 


12 




U2-U3 


26.4 


53 2.63 


14 




#2 night: 2004-07-08 


O dlilL)id.LUiS 






U2-U3 


44.9 


11 1.38 


19 




U2-U3 


43.5 


50 1.11 


14 




#3 night: 2005-05-25 


7 calibrators 






U3-U4 


45.8 


89 1.36 


14 




#4 night: 2005-06-23 


4 calibrators 






U3-U4 


45 


89 1.39 


23 




U3-U4 


46.8 


90 1.33 


16 




#5 night: 2005-06-27 


4 calibrators 






U3-U4 


54 


139 1.29 


11 




#6 night: 2005-07-20 


7 calibrators 






U3-U4 


58.6 


101 1.07 


22 




U3-U4 


59.6 


124 1.08 


22 




#7 night: 2005-08-23 


2 calibrators 






U3-U4 


62.5 


112 1.00 


20 




U3-U4 


54.9 


136 1.25 


21 




Table 2. Parameters of the calibrators used for estimating th 


transfer function. 








Name 


Spectral 


Diameter Flux dens. 


Night " 




type 


[mas] [ly @ 12yum] 






HD 107446 


K3.5III 


4.42±0.03 22.0 




4 


HD109379 


G5II 


3.25±0.02 15.1 




1 


HD123139 


KOIIIb 


5.34±0.03 36.4 




3 


HD 134505 


G8III 


2.50 ±0.01 8.4 




4 


HD 142804 


Ml III 


2.79±0.04 8.0 




7 


HD152820 


K5III 


2.62±0.01 7.5 




3,4 


HD 160668 


K5III 


2.28±0.01 5.7 




3,5 


HD165135 


KOIII 


3.48±0.02 15.5 




1-7 


HD 169767 


G9III 


2.16±0.01 5.9 




3 


HD169916 


Kllllb 


4.00±0.03 20.1 




4 


HD 173484 


M4III 


3.43±0.04 11.9 




5,6 


HD177716 


Kllllb 


3.76±0.03 17.1 




2 


HD178345 


KOII 


2.49±0.01 7.6 




5-7 


HD188512 


G8IV 


2.07+0.01 5.9 




1-3 


HD 192947 


G6/G8III 


2.33±0.02 7.5 




1 


HD2 13009 


G7III 


2.08±0.02 5.9 




5 


HD2 18594 


Kim 


3.18±0.02 12.4 




3 


HD220704 


K4III 


3.45±0.02 13.2 
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The night in which the calibrator was used 



8.2m unit telescope s of the ESO Paranal Observatory in Chile 
(iLeinert et al.ll2003h . We used the standard 0.5"x2" sHt and dis- 
persed the light over the entire A^-band (8-13 /im) with the prism 
providing a spectral resolution of 7? ~ 30. The first MIR fringes 
of IRS 3 were recorded successfully on the night of 7 July 2004. 
The whole dataset comprises 14 independent visibility measure- 
ments of IRS 3 with interlaced calibration measurements. The 
projected baseline length (PB) and position angle (PA) of each 
measurement are given in Table[T]for the beginning of the fringe 
measurement. The last column gives the median relative uncer- 




small errorbars; 
statistical error estimation 
without adequate lowflux 
background uncertainty 



Fig. 3. The influence of faulty background estimation on the fi- 
nal accuracy for data from 2005-05-25. The standard reduction is 
plotted (solid line) with respective uncertainties (small errorbars 
to the right of larger ones). It suggests a significant deviation 
from a smooth, flat visibility spectrum (dotted line: weighted 
quadratic fit over the silicate absorption feature). In particular 
the visibility increase (solid line) at /} < 11 //m towards the ab- 
sorption center cannot be confirmed if we consider the impact of 
background subtraction errors (large error bars). To demonstrate 
the flux dependence of the importance of background accuracy, 
we overplotted the total flux spectrum (in Jy). 



tainty obtained for the caUbrated visibiUty outside the low-flux 
region of the silicate absorption. 

2. 7 . Interferometric calibrator stars 

In Table |2] we list the used calibrators with their main fea- 
tures. The angular diameter ascertainments result from fitting 
Atlas9 and Marcs model atmosphere spectral energy distribu- 
tion s (SEDs) to optical and NIR photome try. The chosen mod- 
els (Kuruczl ll992l 11994; Pl ez et alJlT992l) adopt solar metaflic- 
ity. Details of the model fitting are given in van Boekel (20041 
chapter 5). The modeled llfim flux densit ies listed in Table 2] 
are consistent within <5% with the work of ICohen et al.l (1 19991) 7 
who present a list of absolutely calibrated infrared spectra. 
Furthermore, none of those calibrators shows an MIR-excess, 
defined as having a measured 12 jjm flux density more than 3 cr 
above the synthetic spectra fitted to the optical-NIR data. Such 
an MIR-excess would indicate the existence of (extended) dust 
shells. Dust shells, which can be expected to exist around K-M 
giants, radiate the stellar flux at longer wavelengths and decrease 
the visibility amplitudes of the entire object; i.e., they deteri- 
orate the calibrator properties. The absence of an MIR-excess 
also makes all visibility calibrators usable as potential photo- 
metric calibrators, down to the 5% uncertainty, which is often 
not reached due to atmospheric and instrumental variability. The 
photometric variability and the fitting uncertainty of the angular 
diameter of all the calibration stars in Table|2]affect the accuracy 
of the derived transfer functions by less then 1% at all VLTI 
baselines (< 200 m). 



2.2. Calibration and absolute accuracy 

The immediate measurand, the raw visibility Vi-aw, is the ratio of 
measured correlated (Fc) and total flux (Ft). The final quantity 
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Fig. 4. Transfer functions in gray, derived from the individual 
calibrator measurements of 2005-05-25. The error bars around 
5% indicate the statistical error of the single data frames of one 
measurement, which is slightly smaller than the absolute vari- 
ation in the TF over the night. The mean TF, which we used 
for this night, is overplotted in bold solid black style with the 
standard deviation. The two remote TF (gray), highlighted with 
a diamond symbol, were rejected from calculating the mean be- 
cause of abnormal behavior relative to the average of the good 
TF. 
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Fig. 5. Different analyses of the TF, derived from calibrator 
measurements on 2005-05-25. To increase the clarity of the plot, 
the wavelengths of two overplots have been shifted, as indi- 
cated by the respectively shifted labeling. Left error bars: av- 
erage (mean & standard deviation) of the individual TF. The 
two other data sets are averaged in normalized space (division of 
the individual TF by their integral). Central: maximum differen- 
tial fluctuation between adjacent spectral channels, presenting a 
conservative maximum of the spectral shape error of TF. Right: 
standard deviation of the normalization factor, representing the 
order of magnitude of the variability of transmission between the 
Fc and Fj measurement. 



of interest is the calibrated visibility (Vcai), which is computed by 
dividing the raw visibility (Vraw) by the interferometric transfer 
function (TF) of the observation. Propagation of errors relates 
the uncertainties A to each other: 



calib 



TF = 



^raw 



calib 



ideal 



calibrator 



Vcalib 



\ TF 













/AKawV 







IATFV 



(1) 



(2) 



calibrator 



The TF is derived from the Viaw measurement of a calibrator star 
and its intrinsic visibility (Eq.|2]i. The uncertainty of the diameter 
of the calibrator is usually too small to affect ATF (Table |2]i. 
Thus AVcaiib suffers twice from the accuracy of the estimation of 
Vi-aw (Eq. [TJ. This central accuracy can be estimated by means 
of the following considerations and investigations of calibrator 
data sets. 

All constant defects of Vi-nv/ are corrected for by multiplica- 
tion with TF and thus have only a minor influence on the final 
accuracy; e.g., the time delay between the measurement of Fc 
and Fj is only used a few minutes in the high-sensitivity mode 
of MIDI. This reduces the variation of atmospheric transmission 
in the thermal infrared, due to airmass difference and temporal 
fluctuations, to such a level that a small and time-constant im- 
pact on the quotient Fc/F-y can be assumed. We prove this by the 
following test, visualized in Fig.|5] If the atmospheric and instru- 
mental flux transmission change randomly at a significant level 
between the Fc and Fj measurements, the orders of magnitude 
of Vraw and the resulting TF are affected and the absolute accu- 
racy of the calibrated visibility is worse than the relative spectral 
accuracy. To derive this spectral accuracy we linearly normalize 
each individual TF of a night by dividing through the integral 
over a fixed part of the A^-band and estimate the maximum scat- 
ter in each spectral channel. If this scatter is significantly smaller 
than the scatter of the normalization factors, a strong transmis- 
sion variability has been observed (assuming constant intrinsic 
visibility). Figure |5] shows that this is not the case. Note in the 
same figure the increase in the shape error around the atmo- 
spheric ozone absorption band at 9.7 jum and at the borders of 
the A^-band, where the applicability of simple linear normaliza- 
tion decreases. 

Faulty and unequal background suppression during reduc- 
tion of both data sets also does not vary on a significant level for 
the bright calibrator stars. The influence of background noise for 
faint science targets is di scussed at the end of the section. 

As a matter of fact, iLeinert et al.l (12004,) state that the ac- 
curacy of Vraw is dominated by the accuracy of the overlap be- 
tween the interfering beams, which critically relies on the wave- 
front corrections during observations. This is difficult to quantify 
and thus cannot be easily incorporated into the TF. Usually the 
quality of the beam overlap does not vary statistically during 
one observation. This means that there is no mean overlap accu- 
racy to correctly describe a single observation and that the TF, 
if derived from one calibrator measurement, may not be appli- 
cable to a subsequent science observation. We confirm this on 
the basis of our dataset. A single observation of both Fc and Fj 
consists of several frames that are averaged during data reduc- 
tion. We estimate the standard deviation of the subsequent scans 
and analyze the background in chopped photometry data to de- 
rive a statistical uncertainty for each measurement of Vraw The 
resulting total statistical uncertainty cr(yi-aw) is similar for the in- 
dividual calibrator observations, but it cannot fully explain the 
larger scatter of individually estimated TF over the respective 
observing night. Figure |4] demonstrates that the standard devia- 
tion over several calibrators is only as small as expected from the 
estimated uncertainty of the individual measurements (~ 5%), if 
outliers are rejected from the average. That means that the abso- 
lute accuracy of the visibility measurement can change signifi- 
cantly with each new pointing. 

Furthermore, the overnight scatter of TF is significantly re- 
duced in the 2005 data due to the more stable VLTI feeding by 
the higher-order AO-system Macao instead of the earlier t ip-tilt- 
only correction of the Strap units dArsenault et al.l 120031) . This 
also confirms the dominating influence of varying beam overlap 
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and flux concentration between different pointings on the visi- 
bility accuracy. 

The above given tests and considerations underhne the origin 
of the accuracy of the Vraw measurement and justify the accuracy 
of a single measurement only being given by the TF statistics of 
several calibrator measurements. We calculated the mean and 
standard deviation of the different TF to quantify TF and ATF. 
Note that the presented estimation of errors even holds for an 
extended calibrator with visibilities significantly lower than one, 
provided that the diameter of the calibrator is known at sufficient 
precision. Although Vraw, calibrator will systematically change over 
the night with changing projected baseline length, this effect is 
annihilated by calculating TF (Eq.|2]l. 

Since the number of calibrators per night is typically around 
8, one TF strongly deviating from the others can influence the 
mean and standard deviation significantly. To avoid this situa- 
tion we did not use the median, but rejected the anomalous TF 
manually (on average one per night, Fig.|4|. This has the advan- 
tage that the mean and standard deviation of the resulting sample 
better represents the spectral shape of the TF reducing the prob- 
ability of artificial spectral features in the calibrated visibility 
spectrum. 

To complete Eq. [1] we have to know the relative uncer- 
tainty of Vraw of the measured target. Qualitatively the origin 
is similar, as discussed above for the calibrators, but the flux of 
the science targets can be significantly smaller than the calibra- 
tor fluxes. Since our calibrators cover a range of 5-35 Jy, we 
searched for any flux dependence of AViaw on the basis of our 
data. We checked for every night whether on average the scatter 
of the TF of the fainter calibrators is greater than the scatter of 
the TF of brighter ones, but did not find any such flux depen- 
dence. Thus we assume a flux-independent relative uncertainty 
of Viaw resulting in the uncertainty of the calibrated visibility: 



AVcalib ^ ^2 . 



Vcalib 



TF 



(3) 



This estimation might still underestimate the relative accuracy 
of the final calibrated visibility of the science target, since IRS 3 
was observed with an off-axis AO-guide star, which decreases 
the AO performance with respect to the on-axis guiding on 
the calibrators. We compared the PSF of the photometric data 
of IRS 3 with the calibrator measurements and specified three 
datasets with calibrator-like PSF. These are highlighted through- 
out the analysis and provide the tightest constraints on our inter- 
pretation of the data (cf. Sect.|2!4l&[3]). 

If the correlated flux drops significantly below 5 Jy and ap- 
proaches the sensitivity limits, an increase in the relative un- 
certainties will probably occur due to the increased influence 
of noise; but outside the deep silicate absorption, the correlated 
spectrum of IRS 3 is above this limit. Furthermore, imperfect 
background subtraction has an increasing influence on the final 
accuracy with decreasing source flux. To account for this, we 
estimated the background level by reducing sky-frames without 
the source. Typically the background randomly varies between 
zero and a value close to F\,„, where Fb„ is the sum of the mean 
and standard deviation of all spectral channels of the background 
frames. This leads to a maximum background induced error in- 
terval of 



Vraw.bg 6 



bg 



Fc 



Fj Fj — -Fx.bg 



(4) 



tion center, which appears to be present after applying the stan- 
dard reduction, cannot be verified after incorporating the back- 
ground uncertainty (Fig. [3]). Indeed, not all datasets show such 
an increase after the standard reduction. 



2.3. Data reduction 

We used the data reduction package Mia-i-Ews provided by the 
MIDI consortium^. This package offers two different methods 
for reducing the data: incoherent averaging of the fringe power 
in each spectral bin over several scans, and coherent averaging 
of the single dispersed scans. In the A^-band, the latter method 
can provide the differential phase information in addition to the 
visibility amplitude, if the atmospheric and instrumental delay 
and dispersion have been removed properly. A more detailed 
description of both methods an d the ir rea l ization in Mia-hEws 
is given by iLeinert et al] (|2004|) and iJaffd (|2004 . respectively. 
Since noise always contributes positively to the power spectrum, 
it is not automatically reduced by averaging incoherently over 
several scans. In contrast, the coherent integration reduces the 
statistical noise of the fringe data by averaging, which makes 
this method favorable at very low correlated flux (below about 
1 Jy) and low SNR. Correlated fluxes down to 0. 1 Jy could be 
estimated by this method (Jaffe, private communication). Since 
our MIDI data is usually well above this limit, we found consis- 
tent results of both reduction algorithms throughout the complete 
dataset (Fig.|2]i. 

The constantly changing baseline projection due to earth ro- 
tation limits the integration times of the Tc measurement. To 
reduce the noise level, which is intrinsically high in the thermal 
infrared, only the pixels with the highest SNR should be con- 
sidered by the reduction algorithm. This is achieved by detec- 
tor aperture masks, which can lead to reduced SNR if the cho- 
sen aperture is too large, too small, or misplaced with respect 
to the incident source photons. This effect is strongly enhanced 
by beam distortions and motions that are only partially corrected 
by the AO system, which holds especially for our data since we 
had to lock the AO on an off-axis guide star 35" away from the 
GC. Therefore we reduced the full dataset four times: with both 
coherent and incoherent fringe averaging and by applying two 
different apertures (or detector masks), a wide standard one with 
fixed location and width and a narrower one that is chosen dy- 
namically for the PSF of every observation to trace the beam 
maximum and to have an optimized width (this is the narrow 
mask in Fig.|2]i. 

After calibrating the data and estimating its accuracy as 
described in Sect. 12.21 the four spectra (ycaiib,i..4 ± Aycaiib,i..4) 
were averaged. To estimate the accuracy conservatively, we se- 
lected a maximum AVcaiib including all four uncertainty intervals 
(Fig. 13. 

Note that the apparent correlation between systematically 
higher visibilities and the use of the narrower mask in this fig- 
ure is not a general systematic effect, but depends on the data 
set. Other data show an apparent correlation between visibil- 
ity and averaging method, independent of the chosen mask. 
Consequently we do not prefer one single reduction method over 
another but use the merged data. Furthermore, the differences are 
significant with respect to the intrinsic uncertainties, resulting in 
clearly increased uncertainties of the merged data. The mean fi- 



which dominates the flanks of the strong silicate absorption to- 
wards IRS 3. A putative increase in visibility towards the absorp- 



^ See the ESO web pages for general information 
and links (http: //w ww, eso. org/instnments/midi / ). 
Currently the data reduction package is provided on: 
http : / /www ■ strw . leidenuniv . nl/~nevec/MIDI/ | 
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nal relative accuracy of each merged scientific dataset is given in 
Table [U 

2.4. Photometric Calibration 

MIDI observations provide the astronomer with spectra of the 
total flux of the target in addition to the visibility modulus. 
We flux-calibrated the spectra on the basis of the regularly ob- 
served interferometric calibrator stars, typically late type gi- 
ants (Sect. 12.1b . We fitted an airmass-dependent system re- 
sponse to the calibrator measurements of each night to flux- 
cahbrate the science d ata (for a more detailed description see 
Ivan Boekel et al.ll2005h . Spectra that were obviously faulty were 
not taken into account. Such a time-independent system re- 
sponse model makes it possible to use all the good calibra- 
tor measurements of the night, which is especially favorable if 
the broad band spectrum of the science target is not known. 
Furthermore, the airmass dependence can minimize the impact 
of larger distances between the calibrators and the target. In 
some nights only one star was closer than 10° to the target be- 
cause of scheduling requirements. On normal nights, incorpo- 
rating airmass-dependence reduced the calibration uncertainties 
by up to 5%. But then the model includes uncertainties due to 
varying atmospheric transmission and the instrumental through- 
put over the night, resulting in final photometric accuracies of 
about 5-10%, which dominate the intrinsic uncertainties of the 
used calibrator spectra. 

The Gaussian detector masks (Sect l2.3l l that are typically ap- 
phed do not affect the visibility calculation but only the photom- 
etry. If the science target is not completely unresolved by the 
single telescope PSF (Fig. [T4l i or if such a weighting mask is 
not well-centered on the brightest pixels, the measured flux is 
decreased. We reduced the photometry separately without apply- 
ing any mask to take care of any such bias. Furthermore, it turned 
out that a lot of datasets on the target have at least one beam of 
significantly lower quality than the other one. Since the calibra- 
tor measurements do not show such a strong beam variation, this 
effect is assumed to result from the use of an off-axis AO guide 
star during the observations of IRS 3, decreasing the accuracy of 
the wave-front correction (Sect l2.3l l. A manual selection of good 
datasets facilitates a final photo-spectrometric accuracy of less 
than 10%, which is decreased towards low fluxes due to the re- 
maining background. The result is shown in the lower spectrum 
in Fig. [Toll. 

3. Results 

In this section we present the direct results of our interferomet- 
ric measurements. First the measured and calibrated data are 
shown. Then in the following sections we discuss how the ob- 
served spectra constrain the underlying brightness distribution. 
This discussion outlines the average properties of IRS 3 in the 
MIR at MIDI resolution. 

3.1. Visibility moduli 

The measured visibility moduli are given in Table [3]. Following 
Sect. 12.21 we did not find within the uncertainties any deviation 
from a smooth visibility slope over the full A^-band. Since no 
correlated fluxes have been measured at the center of the silicate 
absorption, we give the mean visibility and its accuracy for two 
adjacent wavelength intervals in Table [3] The fully calibrated 
visibility spectra are shown in Fig.|6l Note the two outliers at the 



Table 3. Measured mean visibility moduli below and above the 
silicate absorption, centered at 9.8 fim. 



Julian date PB " 


PA121 




Visibility 




[m] [deg E of N] 


8-8.7 


yum 


11.8- 


■13 jim 


#1 night: 2004-07-07 












2453194.7 45.7 


46 


0.17 ± 


0.05 


0.27 


±0.04 


2453194.7 42.3 


52 


0.20 ± 


0.06 


0.28 


±0.04 


2453194.8 37.0 


55 


0.24 ± 


0.07 


0.34 


±0.04 


2453194.9 26.4 


53 


0.47 ± 


0.19 


0.53 


±0.10 


#2 night: 2004-07-08 












2453195.5 44.9 


11 


0.22 ± 


0.06 


0.32 


±0.06 


2453195.7 43.5 


50 


0.22 ± 


0.05 


0.26 


± 0.04 


#3 night: 2005-05-25 












2453516.6 45.8 


89 


0.17 ± 


0.03 


0.25 


±0.04 


#4 night: 2005-06-23 












2453545.5 45 


89 


0.10 ± 


0.05 


0.18 


±0.05 


2453545.5 46.8 


90 


0.11 ± 


0.04 


0.22 


±0.04 


#5 night: 2005-06-27 












2453549.8 54 


139 


0.13 ± 


0.02 


0.22 


±0.03 


#6 night: 2005-07-20 












2453572.5 58.6 


101 


0.10 ± 


0.04 


0.21 


± 0.05 


2453572.7 59.6 


124 


0.08 ± 


0.03 


0.18 


± 0.04 


#7 night: 2005-08-23 












2453606.5 62.5 


112 


0.05 ± 


0.02 


0.15 


± 0.04 


2453606.6 54.9 


136 


0.05 ± 


0.02 


0.12 


±0.04 



" PB and PA stand for projected baseline length and position angle 
and characterize the interferometric resolution at the time of the obser- 
vation. 



bottom of the upper and the middle panels. With respect to their 
projected baseline length, they show visibilities too low to be 
consistent with the other data. This is most probably an artefact 
due to bad beam overlap at those observing times. 

3.2. Probing the circular symmetry in uv-space 

Circular symmetry in uv-space implies circular symmetry of the 
brightness distribution. To probe the variation in the measured 
visibility with changing position angle (PA), we compared our 
data at a fixed uv-radius, where most data were obtained. In 
Fig.|2]we show the uv-coverage of the entire dataset, overplotted 
with a ring of constant uv-radius (5 MA). Note that the definition 
of uv-radius (= projected baselinelength/observing wavelength) 
leads to radial lines in the figure and to the situation where along 
the circle datapoints at different PA might be observed at differ- 
ent wavelengths. The uv-radius gives the angular resolution and 
should show similar visibilities at different PA in the case of cir- 
cular symmetry. In Fig. |8]the mean visibilities at /?uv =5 MA 
of each dataset are plotted. The overplotted horizontal indicates 
that the full dataset still conforms to total circular symmetry. 
Although some data points do not perfectly coincide with circu- 
lar symmetry, it has to be remembered that the shown error bars 
of single measurements may be underestimated, since their es- 
timation relies on averaging several calibrator measurements of 
the observing night, leading to an average uncertainty that may 
be exceeded in individual cases. 

On the other hand, the drawn horizontal shows that a slight 
deviation from circular symmetry is possible. The values around 
PA=( 120+30°) appear to lie on average below the values at 
smaller PA. Possible reasons for such circular asymmetry are 
discussed in Sect. 14.41 

But it has to be mentioned that the data points with the best 
photometric quality (used to fit the horizontal in Fig. [S) have 
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Fig. 6. Calibrated visibility spectra ordered following the inter- 
ferometric resolution to show the dependence on the projected 
baseline length (given to the right). Each curve is plotted in a 
different style. For the sake of clarity we plotted only errorbars 
for the solid curves, and the other curves have uncertainties on a 
similar relative scale. The gray dotted lines show the used linear 
interpolation over the deep silicate absorption, where we do not 
have reliable data. Note the different scaling of the panels. 



RS3: uv— coverage 




Fig. 7. uv-coverage of all our observations. The overplotted ring 
indicates a uv-radius of 5 MA. We have investigated the degree 
of circular symmetry along this annulus (Fig. [8). North is up. 
East to the left. 

RS3; uvradius e(4,9, 5.1) 
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Position angle (Eost-of- North) 

Fig. 8. Measured visibilities at a uv-radius of 5 MA. To guide the 
reader's eye we have overplotted a horizontal representing total 
circular symmetry, which was fitted to the data of best photo- 
metric quality in both beams (three thick error bars). If a linear 
correlation with arbitrary slope is fitted to the full dataset, and a 
negative slope shows a ~ 10% smaller reduced which could 
indicate asymmetry or a wavelength-dependent size (Sect. [T2] l. 



nearly identical visibilities at different PA. And a comparison 
with the uv-coverage (Fig.|7]l shows that the longer wavelength 
spectral channels have been considered at those PA with a ten- 
dency towards lower visibilities (Fig. |8j, i.e. the possibly indi- 
cated deviation from circular symmetry in Fig. |8] might in fact 
derive from the slightly wavelength-dependent size of IRS 3. 
We do not expect to observe wavelength-dependent sizes due 
to line emission or absorption of a certain layer, since the vis- 
ibility spectrum over the A^-band is apparently free of spectral 
Hue features within the uncertainties. 

But a slight change in the shape and the size scale of the 
brightness distribution of the dust can be expected due to typi- 
cally lower temperatures of larger, outer dust layers dominating 
at longer wavelengths. This interpretation can easily explain the 
deviations from perfect circular symmetry in Fig.[8]and is further 
backed up by a slight increase of size with wavelength indicated 
by the wavelength dependent analysis of the data (Sect. l3A2l i. 

Circular symmetry implies a vanishing differential phase, 
which remains in the data after coherent averaging and calibra- 



tion. We did not find differential phases along the A^-band spec- 
trum larger than the remaining scatter of less than + 5° around 
zero. Based on the circular symmetry found, we show the PA- 
averaged data plotted over the uv-radius in the lower panel of 
Fig- El Also the black solid line in both panels is the wavelength 
average of all data and lies smoothly between data sets at 8.5 fi 
and 12.5 fim (upper panel). This instead indicates a slight size 
increase in the brightness distribution with wavelength than a 
qualitative change because the overall trend of the visibility with 
baseline length coding the brightness distribution is similar over 
the A^-band. 

Naturally the merged data represent an average bright- 
ness distribution at the cost of losing the slight wavelength- 
dependence but gaining additional and more accurate datapoints, 
since we can use all spectral channels together Note in the lower 
panel in Fig. |9] that statistical data averaging at one uv-radius 
normally reduces the errorbars with respect to the original data. 
The wavelength-averaged data is used only in Sect. l3.4.T] to de- 
rive an overall description of the underlying brightness distri- 
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RS3: wavelength dependent visibilities 



errorbors: 
gray: 8.5^m 
black: 12,5/xm 

black solid line: 
wavelength average 




RS3: Fitted wltln tau(9.8) = 7.16 and T = 410 K 
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uv — radius [MX] 
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Fig. 9. Azimuthal data average in both panels. Upper panel: er- 
ror bars represent the data taken at 8.5 (gray) and 12.5 yum (black 
& diamonds), respectively. In the lower panel, the same spec- 
tral average (black solid line and error bars) is shown with all 
measured visibilities (error bars in light gray). The three datasets 
with the best beam quality are highlighted in dark gray. To cal- 
culate the average, we interpolated linearly over the silicate ab- 
sorption where necessary. 



bution, which itself is used as a basis for describing the wave- 
length dependence of IRS 3 by fitting the average model to the 
wavelength-dependent data (Sect. 13.4.2b . The three datasets of 
highest photometric quality are highlighted. This suggests that, 
beyond ~ 6 MA, the averaged uv-data (black solid line) indicate 
visibilities, which are too low since they lie below the best data. 

3.3. Photo-spectrometry 

In Fig.[10l the flux-calibrated spectrum of IRS 3 is shown based 
on the reduction described in Sect. 12.41 Only the data of best pho- 
tometric quality and AO correction have been taken into account. 
The three resulting data sets, observed in July 2004, May 2005, 
and June 2005, do not show significant photometric variability 
beyond the general uncertainties. This agrees with the study of 
flux-variable sources in the GC bv lOtt et all (Il999h . who find no 
variability for IRS 3. 

3.4. Morphological interpretation of the visibility data 

In this section we explain the MIDI data by a model for the 
brightness distribution, which is as simple as possible, but as 
complex as needed. Although the results cannot compete with 




Fig. 10. Flux-calibrated and dereddened photometry. The up- 
per spectrum is dereddened with rg.g - 7.2, and its error 
bars indicate the wavelength intervals used to fit the temper- 
ature. The middle spectrum is dereddened with Tgs - 3.3, 
which corresponds to the standard average optical extinction of 
A v=25 towards the ce ntral parsec, assuming the extinction law 
bv lMoneti et al.l ( 1200 ih . The lower spectrum is the extinguished, 
measured spectrum. The black solid lines show the extinguished 
and dereddened ;if^-minimized temperature fit of r=410 K. 
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Fig. 11. Best-fit models of the brightness distribution with a sin- 
gle component model; left panel shows a uniform disc, right 
panel shows a Gaussian. The angular diameter is indicated. 



detailed radiative transfer models, they summarize the average 
order-of-magnitude properties in terms of morphological shape, 
size, and flux. Such a heuristic model is therefore an impor- 
tant check and starting point for further analysis, discussion, and 
interpretation of the data (Sect. |4]i. Thanks to the large multi- 
wavelength data set, we get a non-trivial model of the dust dis- 
tribution around IRS 3. 



3.4.1. General structure 

Following the results of Sect. 13.21 we assume a circularly 
symmetric brightness distribution, at first of a wavelength- 
independent size, to derive the general shape of the bright- 
ness distribution. The full dataset is shown in a radial uv-plot 
in Fig. |9] The two simplest, but often applicable, circularly 
symmetric brightness distributions are a uniform disc and a 
Gaussian. Figure [TT]demonstrates that the visibility moduli can- 
not be described by such a simple single component. Despite the 
conservative error estimates, neither model can reproduce the 
data even approximately. 
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Fig. 12. Best-fit models of the brightness distribution with a two- 
component model. Up: two uniform disc components; Bottom: 
two Gaussian components. The error bars are the azimuthally 
averaged data (Fig.|9] Sect. l3.4.TT i. In addition, the probed spatial 
scales are indicated by vertical lines. To the right of each of these 
lines, the visibility of a Gaussian component of the indicated 
FWHM would contribute less than 10% of its flux. The dashed 
line stems from a model adding a third Gaussian of arbitrary size 
larger than 80 mas and illustrates that up to 40% of the total flux 
could have been resolved out by the interferometer. 



In contrast, two superposed components of different sizes 
and flux ratios are sufficient to model the data (Fig.[T2]l. Both two 
uniform discs and two Gaussians fit the error bars satisfactorily, 
but the Gaussians are closer to the measured data: the reduced 
best-fit is about ten times smaller for the two-Gaussian model. 
Furthermore, the data reduction indicates that the visibility at 
low uv-radii (~ 2.5 M.A) is probably overestimated and should be 
expected to lie in the lower half of the indicated error bars. This 
would favor a Gaussian shape for the larger component. But as 
shown by the dashed line in the lower panel of Fig. [12] the inter- 
mediate 50 mas spatial scales are relatively loosely constrained 
to contribute between 20% and 70% of the total flux detected by 
MIDI. Up to 40% of the total flux could have been resolved out 
by the interferometer due to the lack of shorter baseline informa- 
tion. Since with the current visibility dataset there are no means 
to further constrain the outer flux contributions, we stay with the 
two component model representing the data and outline where 
further analysis suggest that indeed some source flux could have 



been resolved out (Sect. |43T l. However, the measured correlated 
flux at the higher spatial frequencies (the smaller component of 
the model), which is most important for the further analysis and 
contributes about one third of the total flux, are not significantly 
influenced by this uncertainty. 

At longer uv-radii the smaller component dominates. A 
closer coincidence with the data suggests a Gaussian shape for 
the smaller component, too. But the final distinction between 
Gaussian and disc shape of the smaller component requires ad- 
ditional data at longer baselines. In the case of a uniform disc, 
the visibility would increase again around ~ 10 M.A as indi- 
cated by the overplotted model in Fig [12] (upper panel). We in- 
vestigated the uv-space > 12 M.A with the longest UT-baseline 
(UTl-4: 130m) without fringe detection. This might support a 
Gaussian shape for the smaller component. But it is also possi- 
ble that the visibility increase, indicative of a uniform disc shape 
of the inner component, was too small to be detected with these 
long baselines. However, this non-detection at the long base- 
line supports the primary finding of our source modeling that 
the spatial scales of the smaller component are also resolvedhy 
the VLTI leading to significant constraints on the physical inter- 
pretation of the data (Sect.HJ. 

At this point there is no indication that the two Gaussian 
components used to represent the data must stand for two phys- 
ically distinct entities, such as two dust layers of different radii. 
At the moment they simply appear convenient for describing the 
data, and the larger Gaussian can be seen as representing the 
wings of the observed brightness distribution. In the later sec- 
tions, physical parameters such as flux-calibrated spectra and de- 
rived temperatures support the idea of several physical compo- 
nents, but a more complicatedly shaped single structure cannot 
be completely excluded. In general the situation is similar to an- 
alyzing the uv-data of radio-inter ferometric observations (e.g. of 
a quasar radio jet in iPott et al1 l2005b). where Gaussian com- 
ponents are used to represent mean properties of the observed 
structure, such as size, flux, and location. To avoid confusion we 
simply keep referring to the two Gaussians as being two compo- 
nents. This understanding is further supported by the later anal- 
ysis. 

3.4.2. Wavelength dependence 

The spectroscopically resolved MIDI data allow the investiga- 
tion of wavelength dependence of the observed brightness distri- 
bution. We base this analysis on the two-Gaussian model of the 
previous section, since it fits the wavelength-averaged data set 
perfectly. A model-independent analysis of the observed bright- 
ness distribution by Fourier-transforming the data (based on 
the fundamental van Citt ert-Zernike theorem of interferometry, 
e.g. iLabevrie et al]|2006l) confirms of our model approach, but 
the limited spatial frequency coverage hampers further model- 
independent conclusions from such an analysis. We describe the 
existing wavelength-dependence of the data in terms of chang- 
ing sizes and relative flux contributions of the two Gaussians, 
keeping the comment at the end of the previous section in 
mind. Furthermore, the change of the fit parameters with wave- 
length turns out to be reasonably smooth (Fig.[T3Tl. This finding 
backs the application of the wavelength-independent model of 
Sect. l3.4.Tl as a basis for fitting the wavelength dependence. 

We fit intensity and size of two superposed Gaussians to the 
data (Fig.[T3]l. We bin the data to a 0.5 //m sampling, using error- 
weighted visibility averages. The unreliable data around the cen- 
ter of the silicate absorption are interpolated using a x'-f^^ of a 
quadratic curve to the spectral channels of good SNR. In Fig. [6] 
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Fig. 13. Wavelength-dependent model consisting of two 
Gaussians. The trends are discussed in the text. Around 10 fim, 
the Gaussians are fitted to the interpolated data points. I^/i de- 
note the intensities (in arbitrary units) and 0s/i the FWHM of the 
smaller and the larger components, respectively. The flux ratio 
calculates as (Is&^)/{Ii6^). In Sect. 14.51 the physical properties 
of the inner component are analyzed by RT calculations. The 
Gaussian parameters have been fitted to the directly measured 
visibilities in the spatial frequency domain. 




Fig. 15. Zoom of the lower spectrum in Fig. [10] into the wave- 
length interval of the 11.3 jum feature of SiC. The measured data 
is plotted in (dashed) gray, and the black solid line represents the 
best-fit extinguished blackbody SED of r = 410 K. 

apparently do not follow the visibility trend of the other curves 
with respect to the baseline lengths). 

The most intriguing result is that the smaller component (0j) 
shows a roughly constant FWHM of about 18 mas, while the 
larger shows a significant linear size increase with wavelength. 
This might be an indication that we directly resolve the inner 
zone of dust formation at all wavelengths, whic h has a fixed size, 
and th at the dust shell might be carbon-rich jlvezic & Elitzurl 
Il996bh . 

Furthermore, the relative flux contribution of the larger com- 
ponent increases with wavelength, suggesting that the larger 
component represents the outer, cooler dust around the central 
object. Again this is a more qualitative statement. The second 
component stands for the more extended flux of IRS 3, its prop- 
erties average over the conditions at the outer parts, but a smooth 
transition between the inner and outer regions around the star 
cannot be excluded. Temperatures and luminosities of dust on 
the size scales of both components are derived in Sect. l4.3l based 
purely on the interferometric data. 



the individual regions of reliable data are indicated by error bars. 
The interpolated data are used in the central wavelength interval 
only, which is devoid of error bars. The respective dependence 
of the visibilities on the position angle is more visible in Fig.|8] 
In Fig.[T3]the resulting best-fit parameters are shown. Using 
;^^-minimization, the uncertainties of the fit are nearly the same 
as the scatter of the data around the overplotted linear correla- 
tions. A slightly larger systematic uncertainty might be intro- 
duced by the estimated errors of the individual data points. To 
address this, we re-fit the data with increased weighting of the 
three datasets of best photometric quality. The overall trends 
are similar, but the size of the larger component in the two- 
component model approach may be underestimated in Fig.[T3]bv 
5-10 mas, which indicates more flux on low spatial frequencies. 
The FWHM of the smaller component, 6s, may show a slight 
size-increase with wavelength by about 3 mas, and the flux ratio 
increases towards the smaller component {Fs/Fi ~ 0.7). This in- 
creased flux ratio can be understood by lower photometric qual- 
ity, typically decreasing the average visibility due to an imper- 
fect beam overlap. But a decreased visibility means increased 
relative brightness of the larger component. Thus, probably a 
few of the data sets of lower photometric quality systematically 
show visibilities that are too small, thus artificially implying 
structures that are less concentrated than the real ones (see the 
bottom curves of the upper and central panel of Fig. |6] which 



4. Discussion 

4.1. Interstellar absorption and the composition of the 
circumstellar absorbing dust 

Following the most recent published results, we assume a spec- 
tral profile o f the interstel l ar abs orption towards the GC, as 
published by iMoneti et al.l (1200 lb . and a n average visual ex- 
tinction of Av ^ 25 towards the GC dScoville et alJ llool 
IViehmann et al.1 12005|) . iMoneti et all (1200 ll) incorporated into 
their model that the mean interstellar dust towards sources in 
the central 2 pc of the Galaxy shows a relativ ely stronger sili- 
cate a bsorption than in the solar neighborhood dRoche & AitkenI 
fT985h . 

In the NIR no strong local increase in the interstellar extinc- 
tion towards t he region of IRS 3 i s found (e.g. most recently 
confirmed bv ISchodel et al.l [2007al) . However, several authors 
claim additional silicate absorption in the A^-band only along 
the line of sight to IRS 3, possibly intrinsic to the source (e.g. 
Roche & Aitken 1985; Viehmann et al. 2006). 

A first glimpse of the probable location of this additional 
extinction is given by our high-resolution single-telescope A^- 
band imaging with the new VLT/VISIR instrument. In Fig. [14] 
the complete emission of IRS 3 (left panel), unresolved by ear- 
lier imaging, is clearly resolved into a diffuse and a compact 
component (central and right panels, respectively). Although the 
total diffuse flux is even larger than the compact flux, its surface 
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Fig. 14. Our VISIR 8.6 yt/m imaging data of a 4" field of view centered on IRS 3. Left: A zoom of Fig.[T] Middle: The diffuse 
emission w ith subtracted PSF. Right: The flux-calibrated PSF, i.e. the unresolved emission seen by MIDI. The data reduction is 
detailed in jSchodel et alJ l20Q7b). All images have the same color coding as indicated at the top and the same contour lines for 
ease of comparison. The logarithmic contours levels are 1.6"-7 mJy. The integrated flux of the different components is given in the 
figures. Note that the diffuse flux peaks below 40 mJy/pix. 



brightness is very low and hidden in the noise of all MIDI data 
due to the shorter integration times. Thus most of the diffuse flux 
is not included in the single-telescope MIDI photometry {Fi x), 
although it is observed at similar spatial resolution to the VISIR 
data. Therefore the flux-calibrated MIDI Fi spectrum (Fig. [TOl l 
fits the unresolved flux shown in the right panel. The two com- 
ponents, fitted to the interferometric data (including the possible 
fraction of fully resolved flux; Sect. |3.4| i and discussed in the 
following sections, together make up the unresolved emission 
in the right panel of Fig. [14] since the single-telescope observa- 
tions provide an angular resolution of about 250 mas. To clarify 
the situation, we speak of local, interstellar silicate absorption 
in addition to the GC average, if the absorbing silicate is located 
in the dust, which radiates the diffuse emission that is directly 
visible only in the VISIR data. In contrast, absorption in the in- 
ner dust components, resolved by the MIDI flux and visibility 
estimates, is labeled as intrinsic, circumstellar absorption. 

Our MIDI spectro-photometry confirms the existence of a 
broad 9.8 jjm silicate absorption feature remaining in the data 
after correction for standard GC extinction. We show the mea- 
sured spectrum in Fig. [10] (lower curve). For the dereddening of 
the spectrum, we used the yuCep emission profile of the s ilicate 
feature, as realized in the extinction law bv lMoneti et al.l (1200 ih 
for lines of sight to the GC. Several authors state that this pro- 
file matches both the local ISM absorption and the GC interstel- 
lar silicate abso rption profile best (altho ugh at different relative 
optic al depths; iRoche & Ai tken 1983, 119851; IChiar & Tieieiisl 
I2006h . 

Our IRS 3 spectrum of a spectral resolution of = 30 shows 
this coincidence perfectly. We reddened a single blackbody of 
variable temperature with an absorption spectrum of the nor- 
malized s hape of the silicate feature seen in emission towards 
juCep (see iMoneti et al.|[200ll) and a variable optical depth rg g. 
The shapes of the observed and dereddened spectra coincide 
with the spectra shown by Roche & Aitken (1985) (Fig. [TOl ). 
The best -fitted parameters are a blackbody temperature of 
T = (410 + 30) and T9.8 = (7 ± 0.5). This temperature re- 



sembles other A^-band measurements (e.g. iGezari et al.l I1985L 
found a color-temperature of about 400 K), but is below the 
800 K derived from K- and L-band data (Moultaka et al. 2004|. 
This indicates that a single blackbody may not be appropriate 
for describing the complete NIR-MIR SED. Although the hotter 
component, dominating the K- and L-band, cannot be resolved 
against the cooler outer dust shell by the single telescope spec- 
trum or image at 10 jim with additional SED information, our 
interferometric 10 //m data alone can resolve it (Sect. 143] ). 

Since the center of the absorption feature is hidden in the 
noise of the background subtraction, a certain level of uncer- 
tainty remains in the estimation of Tg.g, but the spectrophotomet- 
ric quality of the wings is good enough to exclude Tg.g < 6.5. 
Assuming Ay - 25, this means Ay/rgg < 4. This is a re- 
markable result, since it doubles the sihcate MIR optical depth 
towards IRS 3 with r espect to the average o f the GC region 
((Av/t9.8)gc ~ 8-10; iRoche & Aifl^cenl 11985'). which itself is 
twice as deep as Ay/rg g in the solar neighborhood. This is 
shown by the middle spectrum in Fig. [10] which is the mea- 
sured spectrum corrected for standard GC values of extinction 
(Ay - 25 and (Ay/rg 8)gc)- The remaining silicate absorption is 
obvious. The aforementioned authors quantify for the first time 
an extra T9 ^ « 0.8 for IRS 3 in addition to already enhanced the 
GC-average. Although the spectral resolution of both datasets is 
comparable, the spatial resolution of the MIDI photometry data^ 
is increased by at least an order of magnitude. 

In addition, at 11.3 yum a significant drop in the data below 
the fit is suggested (Fig.[T5]l, even on the logarithmic scale shown 
in Fig. [To] This further absorption feature, in addition to the 
dominating broad silicate absorption, can be attributed to SiC, 
which peaks around 11.3 jum and has a much narrower spectral 
width than the interstellar silicate feature. If we exclude the data 
from the wavelength region around 11.3 ;um for the;^'^ minimiza- 
tion, the discrepancy between data and fit around the center of 
the SiC feature becomes even stronger, although the fitted tem- 

which is the VLT 10 yum resolution of ~ 250 mas, not the interfero- 
metric resolution 



12 



J.-U. Pott et al.: The enigma of GCIRS 3 



perature and T9 g remain constant in the given interval of uncer- 
tainties. This further supports the existence of a SiC absorption 
feature towards IRS 3, but the sampling of the applied extinction 
law, which does not contain the SiC feature, is with 0.5 fim at 
MIR wavelengths not high enough to properly sample the SiC 
feature. The kink of the blackbody fit at 11 .5 jum is probably 
an artificial feature. Thus a definitive answer regarding the exis- 
tence of absorbing SiC cannot be given. 

Arguments against an intrinsic silicate absorption in the im- 
mediate circumstellar dusty environment of IRS 3, which could 
be evoked by a deep O-rich dust shell, are: 

- Most silicate-rich dust shells show the silicate feature in 
emission. Similarly the SiC feature is typically found in 
emission in the dust shells of evolved stars. Since our esti- 
mated optical depth rg g is already very deep, an even larger 
amount of absorbing dust would be necessary to overcome 
the circumstellar emission and result in such strong features 
as observed. 

- The spectral shape of the observed silicate absorption coin- 
cides perfectly with the interstellar absorption features. No 
indications of circumstellar crystalline silicates are found, 
although spectroscopic data with higher SNR and spectral 
resolution covering the full A^-band are needed to investi- 
gate the spectral shape in more detail. 

Thus a circumstellar dust shell free of a significant amount 
of silicates appears to be a reasonable assumption for the imme- 
diate environment of IRS 3. This is confirmed by our visibility 
data, which do not show any spectral feature coinciding with the 
broad shape of the 9.8-silicate feature at any baseline length. 

Such a lack of an intrinsic, circumstellar silicate-rich dust 
shell and the deep interstellar silicate absorption would favor the 
bright IRS 3 to be the primary target for estimating the true spec- 
tral shape of the interstellar absorption in the A^-band towards 
central GC sources at the high spatial and spectral resolution 
now available at ground-based 8 m class telescopes. 

4.2. Dust temperatures from the spectral energy distribution 

The MIR regime is dominated by thermal dust emission. We 
have shown in the previous section that the A^-band spectrum 
of IRS 3 can be fitted convincingly by a single reddened black- 
body. But published studies of stellar dust shells show that the 
full infrared photometric information is required to describe the 
optical and physical properties of the shells. Accordingly, we 
investigated the complete infrared wavelength range from 1.6- 
20 /vm, available at a spatial resolution sufficient for distinguish- 
ing IRS 3 from other sources. The corresponding dereddened 
SED is shown in Fig. [161 The bolometric flux ratio indicated in 
the figure is better constrained by the MIDI visibility estimates 
(Sect.ig. 

A second star close to IRS 3 has rece ntly been classified 
as a Wolf-Rayet star (Pau mard et al.ll2006l) . It is unresolved in 
the p ublished medium-resolution NACO data dViehmann et alj 
l2005h . but we confirm this secondary star on high-resolution 
NACO images, showing about 30% and less than 10% of the 
IRS 3 H- and /T-band fluxes, respectively. For our SED fit in 
Fig.[T6]we used the accordingly reduced published H magnitude. 
Since the companion is located about 120 mas east of IRS 3, its 
MIR flux could also contribute to the MIDI data. But a signifi- 
cant contribution should show up as a binary pattern in the visi- 
bility data (Sect. |4~4l ), which was not observed. Thus we assume 
negligible contamination of the flux of IRS 3 longward of 2 fj.m 
by the WR-star, which is further supported by the bluer NIR 



RS 3: infrared SED 



1 000-00 



0-10 



; 




1 1 1 1 1 


; 








T,= {800±150) K - 






T| = {450± 50) K - 


r 




Bolomet 


ic Fluxratio(l/s)=2-40 












f '■' 


, , , 1 





1 10 100 

[micron] 

Fig. 16. SED temperature fit (solid line) as the superposition 
of two spatially unresolved blackbody SED (dashed lines). The 
dotted line close to the cooler and presumably larger T\ com- 
ponent shows a blackbody SED of T = 410 K, which was 
fitted to the A^-band spectrum alone (Sect. 14. It . The possi- 
ble bolometric flux ratio range is given in the plot. Data are 
taken fromfVie hmann et al. (2005, H-, K-, L-, M-bands) and 
IViehmann et al., (,2006. Q-band) and corrected (see text). The 
A^-band data was taken from our MIDI observations at 8 and 
13 micron, outside the silicate absorption feature. The data was 
extinction-corrected for Ay =25. 



SED of the secondary as measured with NACO and expected for 
a hot WR. 

To account for the uncertainties in the amount of interstel- 
lar silicate absorption towards IRS 3 (Sect. 14. It , we used only 
the MIDI fluxes at 8 and 13 yum outside the 9.8 //m absorp- 
tion feature for the MIR-SED. All data were dereddened with 
the Moneti extinction law scaled to Ay =25. As for the A^-band, 
a broad interstellar silicate absorption feature is located in the 
2-band. At 20 fim, only narrow-band photometry inside the sil- 
icate absorption was available. We dereddened this 2-band data 
by an optical depth of tq ~ 3.5, which is derived from scaling 
the interstellar extinction law to fit our measured enhanced T9 g 
towards IRS 3, exceeding the average GC values. At least the 
plotted SED data for A < 13 yum should be free of any significant 
flux contribution of the diffuse VISIR component, because of its 
low surface brightness (Sect. 14. Il l and presumably cool temper- 
ature. Only the VISIR 2-band data may contain a fraction of 
the diffuse flux; although similar to the MIDI spectra, these Q- 
band data do not reach the low noise level of the 8.6 yum VISIR 
imaging data. 

We successfully fit the extinction-corrected SED with two 
blackbody spectra. The lower temperature, which we attribute to 
the outer component (T/ = (440 + 50) K), coincides with the sin- 
gle blackbody temperature fitted to the A^-band spectrum alone 
(Sect |4.1| l. This single blackbody SED is overplotted in Fig. [16] 
which shows immediately that for A < Sfim additional flux by 
hotter dust is required to explain the measured infrared SED. 
This is consistent with the findings of iMoultaka et alj (|2004|) . 
who fitted a blackbody temperature of 800 K to their 2-4 jum 
data. That hotter dust emission appears at shorter wavelengths 
limits the optical depth of the outer cooler dust, assuming that the 
hotter dust is located inside the cooler dust (see also Sect. 14. 3t . 
The calibration of the 2~band emission, apparently too faint to 
fit the model, was very difficult and the deviations can be fully 
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Table 4. Temperatures from the A^-band MIDI data based on 
size and relative flux estimates of the two-Gaussian model and 
T = 0.5 of the outer dust. 



Properties 


Inner (20 mas) 


Outer (50 mas) 




8 /jm 13 yum 


8 yum 13 yum 


F"[}y] 


35 32 


21..55 25..66 


e>^[mas] 


18 18 


40 55 


Tpmtio [K] 


(610±180) 


(460±100) 


^sphere, 8 [K] 


(830±120) 


(480..630±70) 


Tsphere, 13 [K] 


(1020±210) 


(400..580±90) 



" The individual component fluxes and sizes have uncertainties of 
about 15% and 10%, respectively, within the model. The additional un- 
certainty range of the flux contribution of the larger component due to 
the lack of data at short spatial frequencies (Sect. 13.4.11 ) is indicated 
with dots. 



attributed to calibration errors (Viehmann, priv. communication) 
including the uncertain amount of interstellar silicate extinction 
in that band. 

We can state, summarizing these considerations, that no fur- 
ther blackbody component is needed to model the full near- 
and mid-infrared SED. In particular this excludes any signifi- 
cant stellar contribution to the NIR fluxes of the SED, emerging 
from IRS 3 itself or a second star inside the PSF. This leads 
to two possible interpretations: either the inner and hotter dust 
component is optically thick at NIR wavelengths avoiding any 
direct detection of stellar light or the enshrouded star has a very 
hot continuum emission. The absence of any stellar photospheric 
lines in NIR spectra of IRS 3 (Tanner et al. 2006) supports the 
optical thickness of the circumstellar dust at these wavelengths. 

4.3. Dust temperatures from the interferometric data 

In contrast to the previous section, here we present a derivation 
of temperatures from the spatially resolved MIDI observations 
in MIR as a further step in interpreting this (and similar) inter- 
ferometric data. 

The simplest morphological interpretation of the data con- 
sists of two circular symmetric Gaussian components enshroud- 
ing the same central object (Sect. 13.4b . One-dimensional radia- 
tive transfer calculations confirm bell-shaped brightness distri- 
butions of circumstellar dust shells (Ivezic & Elitzur 1996b). 
That we can distinguish at least two dust components suggests 
the outer one is optically thin in the respective wavelengths range 
and physically separated. That is, the observed total flux consti- 
tutes of the sum of the flux of both components, adjusted by the 
optical depth r. From the observed total fluxes and flux ratios, 
we calculate the component fluxes as 

fL,a - F,,, + Fu = e-^-' + (1 - .-■') Fl, 

= (/;^ + l)Fu (5) 

where F^^^ is the total dereddened flux (Sect.|4]l), R the flux ratio 
between the inner and the outer dust shells as observed, and F' , 

s,l 

are the intrinsic flux densities corrected for the optical depth. All 
measurable quantities in Eq. |5] are supposed to be wavelength- 
dependent. In other sections we present strong indications that 
the absorbing silicate is not located in the inner circumstellar 
dust (Sect. |4!2l & 14.51 ). Nevertheless, here we confine the calcu- 
lation to the edges of the A^-band outside the silicate feature to 
minimize the possible corruption of the results by faulty correc- 
tion for the interstellar extinction. 



In Table |4] the component fluxes and FWHM-sizes for the 
two-Gaussian model of Sect. 13.41 and the temperatures derived 
here are given: Tsphere is the brightness temperature of a spher- 
ical blackbody of radius 6/2 at GC distance emitting the ob- 
served flux, and rpratio the blackbody color temperature derived 
solely from the 8-13 yum color of each component. We corrected 
both component fluxes for an optical depth of the outer shell of 
T = 0.5. Since typically the optical depths at 8 and 13 //m are 
comparable, the color temperature does not change with t. 

Furthermore, RT calculations of spherical dust shells as used 
in Sect. 14.51 indicate that a black body spectrum of rpratio often 
convincingly fits the MIR total spectrum of a dust embedded star, 
particularly well for carbon-rich dust shells. But the respective 
T^sphere of such models, calculated again with a Gaussian bright- 
ness distribution approximation, is often significantly higher 
than rpiatio of the same RT model. This is in line with the find- 
ings presented in Table |4] showing that the temperatures derived 
here only give general trends and temperature ranges without 
imposing tight constraints on the physical temperature. 

We find: 

- A reasonable increase in rpratio of the inner component with 
respect to the outer one. This coincides perfectly with the 
analysis of the infrared SED in the previous section. The 
MIR-interferometric data alone show the inner, hotter dust 
component in contrast to spatially unresolved A^-band pho- 
tometry. But the physical properties of the inner dust are bet- 
ter confined by the complete IR-SED and by more detailed 
radiative transfer calculations (Sect l4.5l l 

- Physically reasonable temperatures can only be obtained 
with an optical depth of ts/b ~ 0.5 due to the outer dust in 
addition to the known interstellar extinction Ay =25. This 
backs the large rg § discussed earlier found towards IRS 3 
but suggests at the same time that this additional silicate ab- 
sorption does not occur in the innermost circumstellar dust. 

- Constraining the optical depth is also important for deriv- 
ing the intrinsic luminosity of ~ 5 ■ 10* Lq only from the 
inner component fluxe s. This luminosity is slightly higher 
than earlier estimates dBecklin et al.lll978l) that were based 
solely on spatially unresolved A^-band photometry. Thus it 
can now be excluded that any additional source outside the 
central 20 mas significantly contributes to the flux and heat- 
ing of the dust around IRS 3, a finding confirmed in the next 
section. 

4.4. Circular symmetry 

In their spe ctroscopic survey of the central cluster, 
iPaumard et akl (12006) classify the star IRS 3E, only 120 mas 
east of IRS 3, as a Wolf-Rayet star of type WC 5/6. Since 
Wolf-Rayet stars of this spectral type are known to be (often 
strong) dust formers, it is possible that IRS 3E is still visible 
in the A^-band. Additionally its spectral classification as a 
carbon-rich WC-star does not conflict with our interpretation of 
the spectral data in terms of a lack of silicate emission. 

The visibility pattern of such a binary system with 120 mas 
separation and east-west orientation is shown in Fig.[T7]for sev- 
eral /?uv The shown pattern has been calculated for two stars of 
equal brightness, and both are individually unresolved by the in- 
terferometer. IRS 3 is resolved by MIDI, which would decrease 
the amplitude of the variation shown in Fig.[T7J but the variation 
per PA is defined by the binary separation only; i.e., the imprint 
of a 120 mas binary system should show several ripples over 
180° rotation, which cannot be confirmed with our data (Fig.[8]l. 
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Fig. 17. Visibility patterns for a 120 mas east- west binary at sev- 
eral uv-radii. The individual components are not resolved. 



Furthermore, we did not observe neither photometric variability 
beyond the calibration uncertainties (Sect. |3.3l l or any non-zero 
differential phases (Sect. [3T2l i. 

The simplest deviation from circular symmetry is an ellip- 
tical morphology or, more generally, a brightness distribution 
of different apparent extension in orthogonal directions. The 
analysis of Fig. [8] allows such an interpretation: Towards PA - 
120 ° the measured visibilities appear to be slightly lower than 
in orthogonal direction, necessitating a larger extension of the 
brightness distribution in this direction. Such a lateral contrac- 
tion could be evoked by the strong stellar wind of a nearby star, 
but more probably this deviation from circular symmetry can be 
attributed to a wavelength dependent size (cf. Sect. l3.2l i. 

4.5. Radiative transfer models 

The Gaussian model analysis of the data demonstrated that the 
interferometer resolves the inner dust around IRS 3. Thus we 
apply a physically self-consistent radiative transfer (RT) model 
here to further investigate the physical and chemical properties 
of the inner dust around IRS 3 addressed before by the smaller 
of the two Gaussians. 

Because of the circular sym metry of the sour ce, we use the 
one-dimensional code Dusty (.Ivezic et al.lll999|) . Since IRS 3 
appears to be very luminous, isolated, and surrounded by a lot 
of dust, it is the most reasonable to assume that IRS 3 is a post 
main-sequence star with strong stellar winds and massive dust 
formation. We followed a heuristic approach and calculated four 
distinct scenarios spanning the space of possible stellar param- 
eters: hot and cold, realized by stellar effective temperatures of 
r.,hot = 2.5 ■ 10^ K and r.,coid = 3-10^ K; C-rich with C/O 
abundance ratios beyond 1, realized by a domination of amor- 
phous carbon grains in the circumstellar dust and 0-rich with a 
dust composition dominated by warm amorphous silicates. 

We applied radi al densit y profiles dominated by radiation 
pressure of the star (ilvezic & Elitzui1ll995i) . While the chemical 
composition and temperatures of the dust strongly influences the 
infrared spectrum, the stellar effective temperature T,, and lumi- 
nosity L» scale the physical size of the system in dependence on 
Tdust, the temperature of the sublimation zone, where dust distri- 
bution starts. For the grain si ze distribution we us ed a power-law 
distribution as published by iMathis et alj (Il977h and an upper 
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Fig. 18. Data equal to the one in Fig. [16] but further corrected 
to match the intrinsic emission of IRS 3 and its surrounding cir- 
cumstellar dust shell. The corrections are discussed in the text, 
and lead to the following intrinsic fluxes: 0.05, 0.5, 13.5, 30, 50, 
46, 7.0 Jy at 1.6, 2.1, 3.8, 4.7, 8, 13, 19 fim. The best RT model 
SED were overplotted, the respective optical depths to are 0.4 
for both C-rich models and 1.6, and 1.3 for the cold and hot star 
oxygen-dominated models, respectively. 



limit on the grain size of a = 0.25 yum, which was successfully 
applied in similar experiments. 

To demonstrate the different constraints imposed by the data 
we present our results in two parts. First the measured SED can 
be related to the optical depth r of the dust and to the temperature 
Tdust at the inner boundary of the dust shell. In the next step the 
remaining ambiguities regarding the stellar temperature and dust 
composition can be scrutinized further by comparing calculated 
visibilities of the RT models to the MIDI data. 



4.5.1 . Modeling the spectral energy distribution 

This part is aggravated by the fact that IRS 3 is observed both 
through a large amount of interstellar extinction towards the GC 
and through the local (partly diffuse) dust of non-negligible op- 
tical depth at MIR wavelengths as shown in earlier sections. 
Some uncertainty derives from the fact that the shape of the in- 
terstellar extinction only at wavelengths shorter than 8 yum is 
well obser ved and shown to be re latively constant throughout 
the Galaxy dlndebetouw et al.l2005h . In addition we do not know 
the exact dust composition and resulting shape of the interstel- 
lar extinction at MIR wavelengths within the central parsec and 
around IRS 3 in particular. Thus the dereddened SED as the ba- 
sis of this section is expected to create uncertainties on the order 
of magnitude of tens of percent in flux. 

This naturally affects the accuracy of the estimation of Tdust 
and r, which both determine the shape of the intrinsic SED. 
Similarly we cannot meaningfully investigate the gradual im- 
portance of secondary ingoing parameters like the density pro- 
file, chemistry (e.g. existence and amount of SiC, crystalline 
olivines), etc. These parameters cannot be fine-tuned unambigu- 
ously by comparison with the dereddened stellar spectrum, so 
fixed standard values are used. 

Since the relative depth of the interstellar silicate absorption 
is unclear, we excluded the data around 9.8 fim for the SED fits 
and used only data at the edges of the A^-band. In addition we 
con-ected the SED for using only 40, 30, and 20 % at 8, 13, 
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and 19 jum of the single telescope flux making allowance for the 
found interferometric flux ratios, thereby eliminating the emis- 
sive contributions from the outer cooler dust. Furthermore, we 
increased the optical depth at /I > 8 /im beyond the Ay = 25 seen 
at NIR wavelengths to account for the extinction caused by the 
outer flux seen by the interferometer (Sect. 1431 ). This is realized 
by upscaling the assumed shape of the interstellar extinction to- 
wards the GC (Moneti et al. 2001) to match at /I > 8 fim an addi- 
tional To ss 0.5 (Table|4]i. This opacity would already explain at 
least half the Tg g excess (Fig.fTOll if the same interstellar GC ex- 
tinction shape would again be used for the outer dust addressed 
by the larger Gaussian in our simple two-component model. It 
is conceivable that, by confining this correction to /i > 8 fim, 
we slightly underestimate the dereddened fluxes at shorter wave- 
lengths, which should also be aff'ected by the outer dust. 

The remainig excess extinction at 9.8 fim could stem from 
either fractionally enhanced silicate dust in the central parsec in 
general or in the outer and (partially diffuse) dust around IRS 3 
only. In the next section we aim at investigating whether this re- 
maining T9 8 excess towards IRS 3 could originate from the inner 
dust at 20 mas scales and is thus produced in an oxygen-rich dust 
shell. Finding the latter would suggest IRS 3 to contribute to its 
possibly silicate-enriched environs. 

We decided to fix Tdust for all models at a typical value for 
subliming dust of 1200 K. Tests show that by varying this Tdust 
by up to 300 K, we still reproduce the earlier dust temperature 
findings and generate model SEDs with deviations from the data 
that can easily be explained by slightly varying the optical depth 
in the outer dust of the stellar surrounding. However, sublimation 
temperatures well below 900 K can be excluded. Furthermore, 
the optical depth needed to create model SED that comply with 
the data is not negligible. This in turn enables us to scrutinize the 
dust composition in the next section since the different optical 
properties of carbon and oxygen-rich dust show stronger impact 
in the case of higher optical depths. 

The probed optical depth range covered two orders of mag- 
nitude (ri3 - 0.1. .10). The best-fit SED models under the given 
assumptions are shown in Fig. [18] The main results can be sum- 
marized as follows: 

- the O-rich SEDs show much poorer coincidence with the 
dereddened data at longer wavelengths; 

- the MIR optical depth of the O-rich models is with tq, n ~ 
1 .5 even outside the silicate feature that significantly larger 
than for C-rich dust (rc.n ~ 0.4); this is the only finding 
that would comply with a scenario in which the Tg g excess 
originates in the innermost dust; 

- the photometry data alone can be fitted convincingly by 
single-shell dust models, and for each dust chemistry both 
hot and cold star models show fits of similar goodness 

4.5.2. Modeling the visibilities 

Similar to the flux correction in Sect. 14. 5. T] we have to alter the 
theoretical visibilities of our RT models here to account for the 
outer flux probed predominantly by the short baselines prior to 
comparing the model to the data. To conform this correction to 
the findings of the previous sections, we add to the respective RT 
model a Gaussian component of twice the RT model flux and of 
FWHM of 35, and 45 mas at 8, and 13 /im, respectively. Again 
these parameters reflect average properties of the dust emission 
at lower spatial frequencies and are variable within the ranges 
implied by the findings of Sect. 13.41 
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Fig. 19. Measured visibilities on the edges of the A^-band outside 
the silicate feature. The data at 13 yum are highlighted by stars. 
The respective visibilities of RT models plus outer Gaussian (see 
text) are overplotted as indicated. Quantitative details of the RT 
models are given in Table |5] 
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Fig. 20. Measured visibilities at 11.5 yum in the wing of the sil- 
icate feature. The best cool star RT models plus Gaussian are 
overplotted, and the dotted lines represent the models at 1 1 .5 fim. 
While the C-rich model at 1 1.5 //m lies between the model visi- 
bilities at 8 and 13 /im, as does the data, the O-rich visibilities at 
11.5 fim are lower than ones on the edge of the A^-band at similar 
baseline lengths. 
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C-rich: solid, 0-rich: dashed 



Table 5. Best parameters of the RT models, shown in Fig. [T8]- 
l22l The dust shell was confined to end at a radius of 10'* 
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Fig. 21. Comparison of the A^-band optical depth of cool star 
models showing the significant change in optical depth along 
the silicate features of O-rich dust. 
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<— East [mas] 

Fig. 22. Mean A^-band surface brightness distribution of our cool 
carbon star model of IRS 3. The inner rim, the dust sublima- 
tion zone, is visible as an annulus and resolved by fitting the 
RT model to our MIDI data. At the bottom the radial profile of 
the respective image is shown, and the image and radial plot are 
equally scaled logarithmically. 



With the given data it remains unclear how much of the outer 
flux is indeed present in a second shell at the 40 mas spatial 
scale and how close this dust emission resembles a Gaussian 
brightness distribution. The apparent existence of such a sec- 
ond component might be an artefact of our relatively simple RT 
models. Similar apparently multi-component dust shells around 
the well-studied cool carbon star IRC-i-10216 could be resolved 
into one complex dust shell by more detailed m odeling (e.g. 
IIvezic&Elitzuij|1996at iMen'shchikov et all l200l . But it was 
demonstrated that our approach reflects the data and the spatial 
information therein. Moreover, the spectral shape of the higher 
spatial frequencies of the innermost dust discussed here is not 
significantly influenced by this outer dust. A variable mass loss 
rate of the central star might be responsible for our not observ- 
ing a simple dust morphology that could be explained by a single 
simple dust shell alone. 

Further it is important to understand that any wavelength- 
dependent interstellar extinction only decreases the flux but not 
the flux-normalized visibility, implying that the following results 
are less affected by the uncertainty in the interstellar MIR extinc- 
tion along the line of sight to IRS 3. 



" r, denotes the stellar effective temperature 

* T^dust is the dust temperature at the sublimation zone, visible as ring 
of radius R\ in Fig. 1221 

Rij, is the size ratio of the stellar and dust sublimation radius R, . 
L,,boi labels the bolometric luminosity of the embedded star. 

100% amorphous Carbon grains l lHanneill 19881) 

^ 100% warm amorphous silicate grains iOssenkopf et al.|[l992h 



The investigation of the four best model scenarios fitted to 
the dereddened SED lead to two important results. First the 
aforementioned ambiguity of T, is probed by our interferomet- 
ric data. The of hot stars with the assumed T^^'^^ tend to be 
too large implying visibilities too low to model the measure- 
ments (Fig. [19] & Table |5]). Increased Tjust could alleviate this 
size problem of the hot star models, and the described uncer- 
tainties in the derivation of the intrinsic fluxes and visibilities of 
IRS 3, as well as the scatter in the data, prohibit a definite answer 
for the stellar temperature. But the comparison of the fitted dust 
sublimation radii (Table |5]) with the simple Gaussian analysis in 
Sect. 13.41 also favors the cool star models since such Gaussian 
fits tend to show FWHM twice as large as the subhmati on di- 
ameter of respective RT models dlvezic &Ehtzuij |l996b). And 
dust sublimation radii of 30-40 AU with stellar luminosities of 
5 ■ lO'* L0 found for the cool star models are common values for 
dust-forming AGB stars. Thus for the first time we can present 
an experimental indication of the exclusion of WR-like hot star 
scenarii for IRS 3. 

The second result of this section rules out siUcates as a domi- 
nating constituent of the circumstellar dust shell. The calculated 
visibilities of the competing O-rich and C-rich dust composi- 
tions are rather similar at 8 and 13 jum, far away from the cen- 
tral silicate absorption, although already here the smaller size 
difference of the C-rich models between both wavelengths are 
closer to the wavelength dependence of the data (Fig.[T9l&l20li. 
However, towards the center of the silicate feature, the apparent 
source size of O-rich dust shells increases significantly due to 
the increase in optical depth (Fig. 12111. Figure l20l shows that this 
effect is still visible at 11.5 //m where we have precise visibility 
measurements. A significant size increase (and visibility drop) at 
the wings of the silicate feature as predicated by all O-rich dust- 
shell models that are close to our data has not been observed. 
Actually the visibility data do not reveal any clear spectral fea- 
ture but instead shows a smooth wavelength dependence within 
the calibration uncertainties (Fig.|6]). 

Summarizing the RT analysis of our data, the most probable 
scenario for the nature of IRS 3 is a cool dust-forming carbon 
star In Fig.|22]we have plotted a mean A^-band brightness dis- 
tribution of our best-fit, C-rich cool star model. The inner rim of 
dust sublimation, resolved by our interferometric experiment, is 
clearly visible at a radius of 5 mas in the brightness map. The 
estimated stellar luminosity (Table |5]l points to a cool carbon 
star on the AGB. Recent studies of ca rbon stars find simila r stel- 
lar and circumstellar dust properties tGroenewegenll 1995b . The 
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reasonable assumption of a stellar wind driven dust density dis- 
tribution, valid for windy post-main sequence stars, enables us to 
calculate the mass loss M, termina l outflow velocity Voo, an d an 
upper limit on the stellar mass M. (IIvezic&Elitzuilll996al) . We 
obtain the following results for the carbon-rich cool star set of 
parameters (cf. Table|5]), which describes the observations best: 



M 



lO^^Moyr"'; 



30 km s 



(6) 



The calculations are based on a gas-to-dust mass ratio of 200 and 
a dust grain bulk density of 3 g cm~ ^. The inherent uncertainties 



are discussed bv^ Ivezic et alj (Il999l) . 

The respective optical depth of the carbon-rich dust shell in 
the NIR is T2^m = 6, confirming that the circumstellar dust is 
optically thick in the NIR, as akeady suggested by the analysis 
of the spatially unresolved SED (Sect. 14.21 ). Thus photospheric 
CO bandheads, typical of red giant NIR spectra, cannot be ob- 
served, which explains the impossibility of spectral classification 
of the central ob ject as described by many authors (most recently 
iTanner et al.ll2006.) . 



5. Conclusions 

We have presented new high-resolution 10 ;um data of the enig- 
matic IRS 3 object in the immediate vicinity of SgrA*. The anal- 
ysis and interpretation of the data contribute to the understand- 
ing of both the nature of the embedded star and the location 
of the unusually strong interstellar silicate absorption towards 
IRS 3. The new high-resolution VISIR imaging data clearly re- 
solve IRS 3 into both a compact and a diffuse emission com- 
ponent for the first time. Most probably the deep 10 ;um silicate 
absorption feature towards IRS 3 takes place partially in the local 
interstellar dust, which is responsible for the substantial diffuse 
MIR emission. 

The interferometric data bear convincing evidence that IRS 3 
is a cool star within a carbon-rich dust-shell without significant 
circumstellar silicate absorption. The clear imprint in the circu- 
larly symmetric interferometric data of the existence of more ex- 
tended circumstellar dust was shown. A one-dimensional radia- 
tive transfer model can successfully explain the complete, avail- 
able near and mid infrared data, spatially unresolved by single 
telescope measurements, if outer dust emission on 40 mas scales 
and beyond and non-negligible additional optical depth is ac- 
counted for. We estimate a radius of 35 AU and a temperature of 
~ 1200 K for the inner dust rim around IRS 3. 

With a bolometric luminosity of 5 -lO** and a stellar tem- 
perature of about 3000 K IRS 3 appears to be a cool luminous 
carbon star, most probably in the helium-core burning phase. 
Such mass-losing stars are important for enriching their envi- 
ronment. Only because of the high angular resolution of the in- 
terferometer could we demonstrate that a large fraction of the 
dust is produced by IRS 3 itself, whereas most other MIR bright 
sources in the central parsec have been shown to heat the sur- 
rounding interstellar d ust that has already formed dTanner et alj 
120051: iPott e"tani2006h . Thus IRS 3 appears to be the prime can- 
didate for observing ongoing dust formation in the immediate 
vicinity of a supermassive black hole. 

In addition to the interstellar silicate absorption, we found 
indications of SiC absorption. For the first time the new gen- 
eration of MIR spectrometers at 8m class telescopes offer the 
possibility of studying in detail the spectral properties of the in- 
terstellar absorption towards GC stars at angular resolution high 
enough to account for source confusion and substantial diffuse 
emission in this outstandingly dense and dusty region. Since a 



C-rich circumstellar dust shell is free of any silicate feature and 
IRS 3 is embedded in an exceptionally large amount of ISM, the 
resulting deep absorption of the smooth background continuum 
of IRS 3 appears well-suited to studying the shape of the inter- 
stellar absorption towards the GC and the extraordinary strength 
of the silicate absorption feature in more detail. The absorption 
profile depends on the chemical composition and is needed to 
accurately correct for the ISM extinction, yielding the intrinsic 
spectra of all GC sources. 

On the basis of our results we propose to grant IRS 3 a cen- 
tral role in such new studies of the interstellar silicate absorp- 
tion towards the GC. Furthermore, a longterm MIR variability 
study is proposed to investigate periodical dust formation around 
IRS 3, suggested by this complex dust morphology. 
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